Astronomy <fc Astrophysics manuscript no. aa0930 


February 2, 2008 


(DOI: will be inserted by hand later) 
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Abstract. The availability of a number of new interferometric measurements of Main Sequence and subgiant 
stars makes it possible to calibrate the surface brightness relations of these stars using exclusively direct angular 
diameter measurements. These empirical laws make it possible to predict the limb darkened angular diameters #ld 
of dwarfs and subgiants using their dereddened Johnson magnitudes, or their effective temperature. The smallest 
intrinsic dispersions of tr < 1 % in 9ld are obtained for the relations based on the K and L magnitudes, for 
instance log <9 LD = 0.0502 {B — L) + 0.5133 — 0.2 L or log 6» LD = 0.0755 (V - K) + 0.5170 - 0.2 K. Our calibrations 
are valid between the spectral types AO and M2 for dwarf stars (with a possible extension to later types when 
using the effective temperature), and between AO and K0 for subgiants. Such relations are particularly useful 
for estimating the angular sizes of calibrators for long-baseline interferometry from readily available broadband 
photometry. 
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1. Introduction 

The surface brightness (hereafter SB) relations link the 
emerging flux per solid angle unit of a light-emitting body 
to its color, or effective temperature. These relations are 
of considerable astrophysical interest, as a well-defined re- 
lation between a particular color index and the surface 
brightness can provide accurate predictions of the stellar 
angular diameters. Such predictions are essential for the 
calibration of long-baseline interferometric observations. 
We propose in the present paper new and accurate cali- 
brations of the SB-color relations of dwarfs and subgiants 
based on direct interferometric measurements of nearby 
members of these two luminosity classes. Our primary 
purpose is to establish reliable relations that can be used 
to predict the angular sizes of calibrator stars for long- 
baseline interferometry. 

After defining the surface brightness relations (Sect. [21, 
we discuss in Sect. H3 the sample of measurements that we 
selected for our calibrations (interferometric and photo- 
metric data). Sect, rj] is dedicated to the calibration of the 
empirical SB relations, relative to the color indices and to 
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the effective temperature, for stars of spectral types AO to 
M2. We also derive inverse relations to estimate the effec- 
tive temperature from broadband photometry and angular 
diameter measurements. As the established relations are 
intended to be used primarily to predict angular diame- 
ters, we discuss in Sect. [3] their associated errors in this 
context. In Sect. we search for a possible instrumental 
bias linked to one of the five interferometric instruments 
represented in our sample. Numerous versions of the SB 
relations have been established in the literature, mostly 
for giants and supergiants, and we discuss them in Sect.[7| 
Main Sequence stars are potentially very good calibrators 
for long-baseline interferometry, and we discuss this par- 
ticular application of our SB relations in Sect. [SJ 



2. Direct and inverse surface brightness relations 

By definition, the bolomctric surface flux / ~ L/D 2 is lin- 
early proportional to T c 4 ff , where L is the bolomctric flux 
of the star, D its bolometric diameter and T c g its effective 
temperature. In consequence, F — log / is a linear func- 
tion of the stellar color indices expressed in magnitudes 
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(logarithmic scale), and SB relations can be fitted using 
(for example) the following expressions: 

F B =a (B- V) + b (1) 

F v = ax (V - K) + h (2) 

F H =a 2 (B- H) Q + b 2 (3) 

where Fx is the surface brightness. When considering a 
perfect blackbody curve, any color can in principle be 
used to obtain the SB. The index designates the dered- 
dened magnitudes, and the Oi and bi coefficients repre- 
sent respectively the slopes and zero points of the different 
versions of the SB relation. The particular expression of 
the SB relation Fv(V — R) is also known as the Barnes- 
Evans (B-E) relation, and is historically the first ver- 
sion to have been calibrated empirically (Barnes, Evans & 
Pearson H)7(i). However, the relatively large intrinsic dis- 
persion of the visible B-E relation has led many authors to 
prefer its infrared counterparts, in particular those based 
on the K band magnitudes (A = 2.0 — 2.4 fim), as in- 
frared wavelengths are less affected by interstellar extinc- 
tion. The surface brightness F\ is given by the following 
expression (Fouque & Gieren I1997f) : 

F x = 4.2207 -0.1m Ao - 0.5 log 6 LD (4) 

where $ld is the limb darkened angular diameter, i.e. the 
angular size of the stellar photosphere. 

The linear expressions of the SB can be inverted easily 
to predict angular diameters, and give linear relations such 
as: 

log0 L D = c 1 (V-K)+d 1 - 0.2V (5) 

for the Fyiy — K) inversion. We have in this example: 

c x = -2a x (6) 

di = 2(4.2207- h) (7) 

In the present paper, we will refer to both the direct and 
inverse relations as "SB relations" . 

3. Selected measurement sample 

3.1. Angular diameters 

Over the past two years, sixteen new angular diame- 
ter measurements of nearby Main Sequence and sub- 
giant stars were obtained with the VLT Interferometer 
(Glindemann et al. l2~0TiUl OTMal I2003b|l equipped with the 
fiber-based beam combiner VINCI (Kervella et al. 2000 
2003a). To complement this sample, we have searched 
the literature, and added to our list the measurements 
related to the stars of luminosity classes IV and V. 
Most of the visible and infrared interferometers are rep- 
resented in our sample, with measurements from the 
Nil (Narrabri Intensity Interferometer, Hanbury Brown, 
Davis & Allen irgp7|l . the Mklll (Shao et al. lP3gg|> . the PTI 



(Palomar Testbed Interferometer; Colavita et al. I1999f) 
and the NPOI (Navy Prototype Optical Interferometer, 
Armstrong et al. I1998J) . Our findings originally included 
a few lunar occultation measurements, but they were re- 
jected as they were related to variable or multiple stars, or 
their precision was not sufficient to give them any weight 
in the fitting process. 

To obtain a consistent sample of limb darkened (LD) 
angular diameters we have retained solely the uniform 
disk (UD) values from the literature. The conversion 
of these model-independent measurements to LD val- 
ues was achieved using the linear LD coefficients u from 
Claret (2000), and the conversion formula from Hanbury 
Brown et al. I|1974a|l . These coefficients are broadband, 
single-parameter approximations of the Kurucz ([1992) 
model atmospheres. They are tabulated for a grid of tem- 
peratures, metallicities and surface gravities and we have 
chosen the closest models to the physical properties of the 
stars. We have considered a uniform microturbulent ve- 
locity of 2km.s~ 1 for all stars. This single source for limb 
darkening corrections ensures the self-consistency of our 
final sample. 

3.2. Photometry 

All the apparent magnitudes that we have retained from 
the literature are expressed in the Johnson system. When 
available, we have preferentially kept the uncertainties 
given by the original authors, otherwise we adopted ar- 
bitrarily a conservative error bar. The U band mag- 
nitudes were obtained from Morel et al. ljT3?S| and 
Mermilliod (1986), and we adopted a ±0.02 error. The 
B, V, R and / bands were obtained from several online 
catalogues available through VIZIER (Ochsenbein, Bauer 
& Marcout EOOOj, and we also adopted a ±0.02 uncer- 
tainty. For the J to L infrared bands, references are not 
so easy to find, as many bright stars are unfortunately ab- 
sent from the recent infrared surveys, like 2MASS (Cutri 
et al.HSII|) or DENIS (Fouque et al. l2lTOjl . We have relied 
on the VIZIER database to obtain the infrared magni- 
tudes of our sample of stars. In some cases, the references 
we used are 30 years old, but many of them have small 
and reliable uncertainties. The original references of the 
measurements are given in the footnotes of Table |3| 

3.3. Data selection 

The SB relations rely on the assumption that stars behave 
like black bodies, i.e. that their colors are mainly governed 
by their effective temperature. A severe deviation from 
this assumption will cause a discrepancy between the ac- 
tual flux per surface unit and the temperature of the star. 

For instance, if there is a second, unresolved star near 
the main object, its additional flux will bias the spectral 
energy distribution. For this reason, we have rejected the 
binary and multiple objects for which separate photome- 
try of the components is not available. 
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The presence of warm material in the circumstcllar 
environment can also create an excess at infrared wave- 
lengths. While this signature is most useful for identifying 
the stars surrounded by protoplanetary disks, it creates a 
bias in the measured color of the star. Some of the stars 
we selected are surrounded by debris disks (e Eri, a PsA, 
rCet, (3 Leo), but the contribution of the circumstellar 
material is negligible, even in the infrared J to L bands. 
The material surrounding these stars is very cold and ra- 
diates mostly in the far infrared domain. We have rejected 
the measurement of (3 Pic from Di Folco et al. ( |in prep.] ), 
due to its large uncertainty (~ 10%) and to the relatively 
high density of its edge-on circumstellar disk that could 
cause significant extinction. 

The fast rotating stars can deviate significantly 
from the black body assumption. As demonstrated by 
the VINCI observations of the nearby Be star a Eri 
(Domiciano de Souza et al. 1200311 . the photosphere of these 
objects can be deformed by their fast rotation. This cre- 
ates differential limb darkening between the pole and the 
equator which appear to have different effective temper- 
atures. This makes it particularly difficult to define the 
true photometric solid angle subtended by these objects. 
In addition, many fast rotating stars go through episodes 
of mass loss, that are likely to create a warm circumstel- 
lar environment. The presence of such hot material around 
the star will create a bias in the flux and color of the star. 
For these reasons, we have chosen to reject the known fast 
rotators (vsini > lOOkm.s -1 ) and the Be stars from our 
list. 

The very low mass stars Proxima (GJ 551, M5.5V) and 
Barnard's star (GJ699, M4V) have been excluded from 
our fitting procedures for three reasons. The first is that 
because of their very low effective temperatures the molec- 
ular absorption bands dominate their spectra and lead to 
a significant discrepancy with the hotter stars. Second, 
these stars are variable, presenting occasional flares that 
make it difficult to estimate their magnitudes. Third, they 
present chromospheric activity that could bias their mag- 
nitudes in the U to V colors. However, we have kept these 
two stars on the SB relation figures for comparison. 

The spectroscopic and eclipsing binaries are less use- 
ful for the estimation of the surface brightness relation, 
as it is in general impossible to measure separately the 
magnitudes of these stars with the required precision. For 
this reason, we have not included in our sample the angu- 
lar diameter measurements obtained by spectroscopic or 
photometric methods. For the interested reader, a rather 
complete compilation of the measurements using these 
techniques can be found in the CADARS catalogue by 
Pasinctti-Fracassini et al. (2001). 

3.4. Final sample 

We report in Tables ^ an d HI the complete set of mea- 
surements that we have considered for our fit. In this 
tabic, the angular diameters #ud (uniform disk) and 



#ld (limb-darkened disk) are expressed in milliarcsec- 
onds (mas). The limb darkening conversion coefficient 
k = #ld /$ud was computed for each star based on the ta- 
bles of Claret l|2000J) . When a physical parameter was not 
available in the literature, it has been estimated approxi- 
mately, and appears in italic characters. The observation 
wavelength A is given as either the name of the photomet- 
ric band (V, H, K) or the actual wavelength in /zm. The 
error bar in the angular diameter of the Sun (G2V) has 
been set arbitrarily to ±0.1%. The parallaxes are from 
the Hipparcos catalogue (Perryman et al. I1997D . except 
the aCen value that was taken from Soderhjelm {1999), 
who derived it from reprocessed Hipparcos data. The in- 
terferometer used for each measurement is indicated in the 
"Instr." column. 

4. Surface brightness relations 

4.1. Fitting procedure 

For each angular diameter measurement #ld, and based 
on the observed apparent magnitudes m\, we have com- 
puted the surface brightness F\ in all bands, using the def- 
inition of Eq. The resulting Fx values were then fitted 
relative to the colors (Co — C\), using a linear model. This 
fit was achieved using a classical \ minimization taking 
into account the errors in both the colors and F\. The 
minimized quantity, using the slope a and zero point b as 
variables, is the reduced \ 2 expression: 



v [(FoY-aiCo-CA-bf 



N -2 



(a Fl ) 2 +a 2 {<j C i) 2 



(8) 



where we have: 



— N the total number of measurements in our sample, 

— (Fo)i the surface brightness of star i in band Co, 

— (Co — G\)i the color of the star of index i computed 
between bands Co and Ci, 

— ere i the 1 a error bar in the chosen color of star i, 

— dpi the lcr error bar in the surface brightness Fq. 

The 1 cr errors a a and a b are subsequently estimated from 
the best fit values a and b by solving numerically the ex- 
pression: 



X 2 rcd (a + a a ,b + a b ) = Xred( a > & ) + L 



(9) 



The solutions of this equation correspond to an elliptic 
contour, due to the correlation between the a and b vari- 
ables. It has to be projected on the a and b axis to give the 



errors. The residuals Afi = Fj — F„ 



Kiel 



are used to esti- 



mate the intrinsic dispersion (7j nt (F) of the surface bright- 
ness relation from: 



1 

N 



N 



J2\(AF t f~(<j Fl ) 



(10) 



This process gives a total number of 72 (a, b) best fit pairs, 
with their associated errors (a a , o~b), and the intrinsic dis- 
persion (7i n t of the data around the best fit model. 
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Table 1. Angular diameters of dwarf stars (luminosity class V) measured by long-baseline interferometry (apart from 
the Sun). They are expressed in mas, and T c g is in K. "Ref.i" designates the reference used for T e g, logg and [Fe/H]. 
When unavailable, the metallicity has been set arbitrarily to the solar value. "Ref.2" designates the reference used for 
each angular diameter measurement (expressed in mas). The errors are given in superscript close to each value. 
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de Strobel et al. I|2001|l : (i) Di Folco et al. {in prep.) ; (j) Morel et al. I|2001|l : (k) Morel et al. (120001 . 

— Ref.2 for angular diameters: (1) Hanbury Brown, Davis & Allen (I1974bt : (2) Kervella et al. Il2004at : (3) Nordgren et 
al. ll2Wrt : (4) Mozurkewich et al. 11217031 : ( 5) Thev enin et al. Jin prep.> ; (6) Bode n et al. (THTOol : (7) Nordgren et al. (tl*9^9l : 
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From there, we can invert these relations easily to ob- 
tain their angular diameter counterparts: 

lo g e LD = c(Co-Ci)+d-0.2C (11) 

The slopes and zero points are computed from the (a, 6) 
pairs through: 

c = -2a, (T c = 2(r a (12) 

d= 2 (4.2207-5), <r d = 2a b (13) 

and the intrinsic dispersions <Ti n t(log Old) are given by: 
^int(log0LD) = 2a int (F A ) (14) 

The same method was used for the fits using the effec- 
tive temperature, except that no error bar was considered 
on the Toff values from the literature (equal weighting), 



and we used a second degree polynomial model instead of 
a linear one. We minimized numerically the following Xrod 
expression using a, b, c as variables: 

Xied {a,b,c) = — — 2^ -g (15) 

JV 6 i=l °F 

where 

-Fmodel(T ff )i = a (logT efr )f + 6(l0gT e ff) i + C. (16) 

The errors in each of the a, b and c coefficients were not 
computed, as the correlations existing between these co- 
efficients make it very difficult to determine them accu- 
rately. This is justified by the fact that the astrophysical 
dispersion of the measurements is largely dominant over 
the 1 <t fitting errors of the model, and the systematic 
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Table 2. Angular diameters of subgiant stars (luminosity class IV) measured by interferometry. The references and 
notations are given in Tabled 
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errors in these coefficients can thus be considered negligi- 
ble. The inversion of the resulting T c g based relations is 
straightforward. With an expression of the form: 



3.00 



lo g LD = d(logT cff ) 2 - 
we have by definition: 
d=-2a, e = -2b 
f = 2 (4.2207 -c) 



e(logT cff ) + /-0.2C A 



(17) 

(18) 
(19) 



As in the previous case based on colors, the intrinsic dis- 
persions CTi n t (log #ld) of the angular diameter relations are 
given by: 



CTint (log Old) = 2 CT in t (F) 



(20) 



In some cases, we could derive only upper limits of 
the intrinsic dispersion CT; n t, as it was found to be smaller 
than the average error of the measurements (in such cases, 
Ea.llUI gives a negative value for CTj n t). For these relations, 
such as log 0LD(2cff , L), we conclude that the intrinsic dis- 
persion is undetectable at our level of sensitivity. 

4.2. Angular diameter relations based on colors 

The SB relations for UBVRIJHKL colors are listed in 
Table H They take the form: 



logf? LD (C , Ci) = c A (C -C 1 )+d x - 0.2 Co 



(21) 



where Co and C\ are any two distinct colors of the Johnson 
system. In many cases, the dependence of the zero magni- 
tude limb darkened angular diameter (ZMLD), defined for 
Co = 0, as a function of the color is not linear in reality. 
Thus, the linear model that we fit does not represent the 
observations well. In this case, we have added a note " nV 
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Fig. 1. Linear fit of the surface brightness relation 
logZMLDs(B — L) (upper part), and the corresponding 
residuals (lower part). The intrinsic dispersion in the re- 
lation is ±0.004 on log ZMLD, equivalent to a systematic 
error of less than 1 % in the predicted angular diameters. 
The open circles designate GJ 699 and Proximo,, which 
were excluded from the fitting procedure. 



after the obtained residual dispersion. The non-linear rela- 
tions should preferably not be used for predictions, though 
the stated dispersions include the non-linearity. 
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Table 3. Apparent magnitudes of the dwarf stars (upper part) and subgiants (lower part) of our sample. The 
uncertainty adopted for each apparent magnitude value is given in superscript. 
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— References: 

(a) Morel et al. (TH?78l . Mermilliod ltT9^6l . 

(6) Morel et al. ltT?78l . Hoffleit & Warren itT^Tl . Perryman et al. itTWt . 

(c) Morel et al. (JTHZHJ, Bessell et al. P^ . Ducati et al. P^u^jl . 

(d) Glass IH97H , Glas s JT9751 . Mould & Hyland (tT5?oTl . Morel et al. iTWM . Leggett lIT^zTl . Ducati et al. tzTTOzTi . Kidger & 
Marti'n-Luis (120031 . the 2MASS catalogue (Cutri et al. 120030 . 

(e) We refered to Colina et al. 119961 for the apparent magnitudes of the Sun. 



In theory, there should be a perfect diagonal symmetry 
between the dispersions listed in Table 0] In reality, the 
symmetry is only approximate, because Co and C\ are not 
symmetric in the expression of log 6*ld(Co, Ci). Therefore, 
an increased dispersion of the apparent magnitudes in one 
band C\ will be reflected preferentially in the #ld(Co, C\) 
dispersion, rather than in 0ld(Ci,Co). For this reason, 
we provide both versions in Table 01 including the quasi- 
symmetric pairs. The best relations based on the K band 
(showing residual dispersions below 1% on the angular 
diameter #ld) are the following: 

log 9 LD = 0.0535 (B — K) + 0.5159 — 0.2 K (22) 
log 6 LD = 0.0755 (V — K) + 0.5170 — 0.2 K. (23) 
and the best relations for the L band are: 
log 6» LD = 0.0412 (U - L) + 0.5167 - 0.2 L (24) 



log 6> LD = 0.0502 (B - L) + 0.5133 — 0.2 L (25) 

log(9 L D = 0.0701 (V-L) + 0.5139- 0.2 L (26) 

log6>LD = 0.1075 (R — L) + 0.5128 — 0.2 L. (27) 

These expressions are valid at least over the range of col- 
ors defined by our sample (Tables ^ an d 121 • I 11 terms of 
spectral types, the angular diameter predictions can be 
considered reliable between AO and M2 for dwarfs, and 
between AO and KO for subgiants. There are indications 
(Fig. that the infrared relations are valid down to the 
spectral type M4V of GJ699, but show some discrepancy 
for the M5.5V star Proxima. The established relations are 
likely to be valid also for subgiants of spectral types later 
than KOIV, but this cannot be verified from our sample. It 
should be stressed that they are applicable only to single 
stars, and the presence of a non-resolved stellar companion 
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Fig. 3. Second degree polynomial fit of ZMLDi(logT c ff) 
(upper part) and the corresponding residuals (lower part). 
The coefficients are given in Table \5\ and correspond 
to a relation of the form logZMLD^ = d(logT c ff) 2 + 
e(logT c ff) + /. The open circles designate GJ699 and 
Proximo, which were excluded from the fit (see Sect. I3.3T) . 
though they are consistent with the model within their 
error bars. 



contributing a significant fraction of the measured flux will 
bias the predicted angular diameters. As more than half 
of the Main Sequence stars are binary or multiple stars, 
care should be taken in the application of these relations. 

4.3. Angular diameter relations based on effective 
temperatures 

Table [5] gives the best fit model coefficients for the rela- 
tions #LD(T e ff, C A ), defined as: 



log LD = d (log T off ) 2 + e (log T cff ) + / - 0.2 C x 



(28) 



The smallest residuals are obtained for the relations based 
on T e ff and the K or L magnitudes, with an upper limit 
on the 1 tr dispersion of 1.0% (the true dispersion is unde- 
tectable from our data): 

log 6> LD = 0.8470 x 2 - 7.0790 x + 15.2731 - 0.2 K 

log 9 LD = 0.6662 x 2 - 5.6609 x + 12.4902 - 0.2 L 

where x = logT c ff. The range of validity of the T e g based 
relations is 3600-10000 K for dwarfs, and 4900-9500 K for 
subgiants. As shown in Fig. there are indications that 



Table 5. SB relations using the magnitude C\ and 
the effective temperature T c g of the star to obtain 
the limb darkened angular diameter #ld (in mas): 
log0 LD (T cff ,C A ) = d(log T eff ) 2 +e log T cff + /-0.2C A . 
The 1 <7 residual dispersions are given in percents of the 
LD angular diameter. 
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the infrared relations are valid for dwarfs with T e ff down 
to -3000 K. 

4.4. T cS (ehD,mx) relations 

By inverting the relations established in Sect. 14.31 it is pos- 
sible to predict the effective temperature of the observed 
stars based on their angular diameter and broadband mag- 
nitude in a single band. As in the previous sections, we 
assume zero interstellar extinction, and the relations are 
valid only for dereddened magnitudes. The formulation of 
the T ff (#LD> m x) laws is easily derived analytically. From 
Eq. [2H1 we obtain log T c ff through the expression: 



V4dlog0LD+O.8dC A 



log T eff = 
that can be rewritten as 

log Teff 

where 
1 



4#-e 



2d 



y/g log 6 LD + hC\ + i + j 



g= d> 



e 

id 2 



0.2 

/l= -' 



/ 
d' 



J 



—e 
2d' 



(29) 



(30) 



(31) 



(32) 



The intrinsic dispersion of the logT ff relations can be 
approximated from the intrinsic dispersion of the log #ld 
relations, as we have (Tint (log #ld) "C 1: 



(Tint (log T e ff) = 0.5 y/g (Tint (log L d) 



(33) 



The corresponding coefficients and dispersions are given in 
Table El The K band relation presents the smallest intrin- 
sic dispersion (er < 0.60%), corresponding to a systematic 
uncertainty of less than 40 K in the predicted temperature 
of a G2V star: 



log Teff = 4.1788 - \/l-18061og(9 L D + 0.2361 K - 0.5695 
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Fig. 2. Johnson B band ZMLD^ relations as a function of color. The errors bars have been omitted for clarity, and 
the fitted models are represented alternatively as solid and dashed lines. From left to right, using the colors: (B — U), 
(B -V), (B - R), (B - I), (B - J), (B - H), (B - K), (B - L). A clear non-linearity is visible on the (B - U) based 
relation. 



2.50 
2.30 
2.10 
1.90 

— 1.70 

to 

CO 

■i- 1.50 
Q 

| 1.30 

o> 
o 

- 1.10 
0.90 
0.70 
0.50 

































■ 














■ 


A 


























■ 




4 x/ 




























a ^ * x, 




+ 












n " * 

* |» * , . V ♦ * m 




□ a 































0.30 



4.00 



3.95 



3.90 



3.85 



3.80 



3.75 
log(Teff [K]) 



3.70 



3.65 



3.60 



3.55 



3.50 



Fig. 4. ZMLDa relations as a function of the effective temperature. The error bars have been omitted for clarity, and 
the fitted models are represented alternatively as solid and dashed lines. From top to bottom, using the zero magnitude 
reference colors U, B, V, R, /, J, H, K and L. 



However, we would like to stress that the uncertainty in 
the measured apparent magnitudes can easily be domi- 
nant, as a ±0.03 error in K will translate into a ±1.7% 
error in T e g, nearly three times as large as the intrinsic 
dispersion. 

Considering that photometry at an absolute level of 
±0.01 is not available for all stars, the T e g predictions from 



different bands can be averaged, taking carefully into ac- 
count the statistical and systematic errors of each relation 
used, in order to reach the intrinsic dispersion level. In ad- 
dition, such an averaging process should not be done for 
stars affected by interstellar or circumstellar extinction, as 
it will affect differently each photometric band. 
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Table 4. Surface brightness relations using UBVRIJHKL based colors to obtain the limb darkened angular diameter 
#ld (in mas) as a function of the magnitude and color of the star through: log#LD(Co, C\) = c\(Cq — C\) + d\ — 0.2 Co- 
The residual dispersions are given in percents of the LD angular diameter. The 1 a errors in each coefficient are given 
in superscript, multiplied by 1000 to reduce the length of each line, i.e. 0.5822 3 2 stands for 0.5822 ± 0.0032. When the 
data depart significantly from our linear fit and present a detectable non-linearity, the dispersion is mentioned in italic 
characters, and we have added the note " nP . The dispersions smaller than 5% are mentioned in bold characters: they 
mark the relations that are the most suitable for predicting stellar angular sizes. 
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C QQO/ ml 
. 00 /0 TLL 


O OA OX 

0.044 /o 


0.0U /o 


Cj 


-0.0818 ' 6 


-0.1043 9 


-0.1581 1 ' 4 


-0.2842 3 2 


-0.7198 169 




0.6280 160 


0.5214 106 


0.4840 9 - 4 


dj 


0.5325 1 9 


0.5216 2,2 


0.5276 24 


0.5299 3 2 


0.5197 6,5 




0.4990 5,8 


0.5066 4,9 


0.5078 4 8 


crj 


5.13% nl 


3.12% 


4.98% 


7.26% nl 


10.9% nl 




10.44% nl 


5.86%, nl 


4.08% 


cu 


-0.0513 ' 5 


-0.0625 ' 6 


-0.0892 09 


-0.1376 1 8 


-0.2290 40 


-0.4312 121 




1.8747 58 - 8 


1.1714 27 ' 7 


du 


0.5189 16 


0.5138 1 ' 8 


0.5145 19 


0.5138 2 3 


0.5093 29 


0.5013 4 ' 4 




0.5352 75 


0.5271 6 - 2 


07/ 


2.67% 


1.24% 


1.12% 


2.06% 


5.53% nl 


10.31% nl 




17.2% nl 


13.0%, nl 


CK 


-0.0440 ' 4 


-0.0535 06 


-0.0755 08 


-0.1144 14 


-0.1805 2 8 


-0.3192 74 


-1.7040 55 1 




3.0857 33 ' 9 


d K 


0.5202 1 - 6 


0.5159 1 ' 6 


0.5170 17 


0.5168 19 


0.5149 2 ' 3 


0.5089 34 


0.5336 70 




0.6258 3 6 


Ok 


2.58% 


< 1.00% 


< 1.00% 


1.60% 


3.67% 


5.87% nl 


1 7. 6%> nl 




26.9%o nl 


cl 


-0.0412 ' 5 


-0.0502 6 


-0.0701 8 


-0.1075 1 4 


-0.1696 2 9 


-0.2826 64 


-0.9843 24 1 


-2.8950 115 7 




di 


0.5167 17 


0.5133 18 


0.5139 19 


0.5128 2 1 


0.5101 25 


0.5081 3 ' 3 


0.5245 5 ' 4 


0.5309 124 






< 1.00% 


< 1.00% 


< 1.00% 


< 1.00% 


2.43% 


3.61% 


13.1% nl 


26.3%> nl 





Table 6. T e g (^ld, m\) relations to obtain the effective 
temperature: logT c ff = — \/g log 6>ld + h C\ + i + j. The 
1 a residual dispersions are given in percents of the effec- 



tive temperature T c g (expressed 


in K). 




C* A 


a 


9 


h 


i 


3 


U 


1.19% 


0.1773 


0.0355 


-0.0703 


4.1186 


B 


1.57% 


0.2721 


0.0544 


-0.0850 


4.2128 


V 


1.61% 


0.3288 


0.0658 


-0.1249 


4.1870 


R 


1.57% 


0.4674 


0.0935 


-0.1848 


4.2120 


1 


1.13% 


1.0155 


0.2031 


-0.2788 


4.4675 


J 


0.60% 


1.0419 


0.2084 


-0.3975 


4.3472 


H 


0.63% 


0.8559 


0.1712 


-0.4230 


4.1149 


K 


< 0.60% 


1.1806 


0.2361 


-0.5695 


4.1788 


L 


< 0.60% 


1.5010 


0.3002 


-0.6977 


4.2486 



4.5. Metallicity 



A possible source of natural dispersion of the SB rela- 
tions is the presence of deep absorption lines in the spec- 
tra of the stars. This effect is stronger for stars that have 
a high metal content. However, as shown in Fig. there 
is no clear evidence of a correlation between the residu- 
als of the least dispersed relation #ld(£, B — L) and the 
metallicity [Fc/H]. This is an indication that our SB rela- 
tions are valid at least for metallicities between —0.5 and 
+0.5 dex, and probably also for lower values. The two 
metal-deficient stars of Fig.0are GJ 15 A ([Fc/H]= -1.40 
dex) and GJ699 ([Fe/H]= —0.90 dex, not included in our 
fits). For typical stars of the solar neighborhood our rela- 
tions are thus always applicable. 
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Fig. 5. Residuals of the fit ofZMLD^L, B — L) as a func- 
tion of the metallicity [Fe/H] of the star. Dwarfs are rep- 
resented by open squares, and subgiants by solid dots. No 
correlation is visible among the stars of our sample. 

5. Sources of uncertainty in the predicted 
diameters 

Several observational and astrophysical sources of uncer- 
tainty add up to create the total error in the predicted 
angular diameters: 

— intrinsic dispersion of the empirical relation: as dis- 
cussed above, the best relations have intrinsic disper- 
sions below 1%. It should be stressed that their predic- 
tions cannot be averaged to reduce this systematic un- 
certainty. However, the predictions from independent 
colors (such as V — K and B — L) can be averaged 
to reduce the statistical uncertainty in the predictions 
due to the errors in the photometric measurements. In 
this process, the systematic uncertainties of each rela- 
tion cannot be reduced and have to be carefully taken 
into account. This is essential as the calibrations have 
been obtained for all colors from the same sample of 
stars, and the resulting systematic errors are therefore 
highly correlated. 

— uncertainties in the apparent magnitudes: combining 
the high precision magnitudes available in the visible 
with the infrared magnitudes produced by the 2MASS 
and DENIS surveys should make it possible to retrieve 
the visible-infrared color indices with a precision bet- 
ter than ±0.02 mag. However, we would like to stress 
that the uncertainty in the apparent magnitude mea- 
surements can easily be the largest contributor to the 
predicted angular diameter errors. The true errors in 
the photometric measurements have to be estimated 
accurately in order to obtain reasonable uncertainties 
in the predicted angular sizes. 

— interstellar extinction and circumstellar matter: for 
the sample of nearby stars that was considered for 
our fits, the interstellar extinction is negligible: apart 



from 7 Gem at 32 pc, all stars are located closer than 
15 pc. However, our SB relations are strictly valid only 
for extinction-corrected magnitudes. The uncertainty 
in the assumed color excess E(B — V) (for instance) 
will translate into an additional uncertainty in the 
dereddened magnitudes. The presence of a significant 
amount of circumstellar matter around the star will 
also affect its spectral properties, and can be difficult 
to detect. 

6. Comparison between interferometers 

The residuals of the fits of the least dispersed relations 
(based on infrared colors) allow us to examine if system- 
atic discrepancies are detectable between the five interfer- 
ometers represented in our sample. For each instrument we 
have computed the average residuals of its measurements, 
and the 1 a error resulting from the averaging of their re- 
spective errors. The results are presented in Table 

We observe that the average residuals are below 1.5 % 
in terms of ZMLD for all instruments. In addition, all 
the deviations are below 1 er, and can therefore be fully 
explained by random statistical dispersion. As a remark, 
the agreement between the VLTI/ VINCI results and the 
Mklll is remarkable, with no systematic deviation de- 
tectable at a level of a few tenths of a percent. This is 
especially encouraging as these two instruments are ob- 
serving at very different wavelengths (visible and K band, 
respectively). 

This comparison exercise relies implicitly on the as- 
sumption that the considered ZMLDo(Co — G\) relations 
are applicable to each instrument's subsample of stars, 
down to the precision of each individual measurement. 
This may not be the case for all stars, but the agreement 
that we observe is a worst case, and the true agreement is 
in any case very satisfactory. 

7. Previous calibrations and other luminosity 
classes 

Previous calibrations of the SB relations for dwarfs 
have been derived by Di Benedetto l|1998|) and Van 
Belle (1999a). These two authors relied on the limited 
sample of hot dwarfs observed with the Narrabri inten- 
sity interferometer (Hanbury Brown et al. 197 lal 1797 lhli . 
The agreement of our calibration with the work by Van 
Belle (1999a) is satisfactory within 1 a for the (V, V — K) 
relation, but there is a difference of about 2 a in the slope 
of the (B, B — K) relation. As the fit obtained by this au- 
thor is based on a small range of colors, we attribute this 
moderate discrepancy to an underestimation of the true 
error bar in the slope, even in the restricted quoted range 
of validity (-0.6 < B- K < +2.0). 

Several calibrations of the SB relations for giants have 
been proposed in recent years, thanks to the availabil- 
ity of a number of direct interferometric measurements 
of this class of stars. Van Belle l|1999a|l used a sample of 
190 giants, complemented by 67 carbon stars and Miras 
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Table 7. Comparison between interferometers. The average residuals of the fits of ZMLDs and ZMLDy for the K 
and L based colors are given together with the corresponding 1 a error bars. All values are expressed in percents of 
the ZMLD values. All the residuals are compatible with zero within their 1 a error bars. 



Instrument 


N 


AZMLD S (B - K) 


AZMLDv(V- K) 


AZMLD B (B - L) 


AZMLDv(V — L) 


PTI 


5 


+0.48 ±0.89% 


±0.82 ±0.94% 


±0.63 ±0.89% 


±0.89 ± 0.93 


Nil 


5 


+1.02 ±1.42% 


±1.03 ±1.46% 


±1.03 ± 1.58% 


±1.02 ± 1.63 


MkHI 


11 


-0.11 + 0.56% 


-0.12 ±0.57% 


-0.06 ±0.64% 


-0.11 ±0.66 


NPOI 


5 


-1.20 + 1.20% 


-1.23 ± 1.24% 


-1.23 ± 1.24% 


-1.39 ± 1.53 


VLTI/ VINCI 


16 


-0.05 ± 0.49 % 


-0.07 ± 0.51 % 


-0.07 ±0.51% 


-0.11 ±0.68 



measured with the PTI (Van Belle et al. [1999b), IOTA 
(e.g. Dyck et al. H998[l and lunar occultation observations 
(e.g. Ridgway et al. H552) l to calibrate the F V (V - K) 
relation of giant and supergiant stars. Welch (1994) and 
Fouque & Gieren (j!997p proposed a calibration of the SB 
relations of Cepheids based on an extrapolation of the 
corresponding relations of giants. Among the supergiants, 
Cepheids occupy a particular place. The observations of 
these variable stars by interferometry, intended primarily 
to study their pulsation, have resulted in the measurement 
of several of these objects (Mourard et al. 1997 Lane et 
al. EST)?]?)! Nordgren et al. IM)T)1 Ke rvella et al. IMiTl Lane 
et al. 120021 Kervella et al. l2004b|> . Based on these obser- 
vations, Nordgren et al. |2002) have established dedicated 
SB relations for Cepheids, and they find a satisfactory 
agreement with previous works. From these studies, it ap- 
pears that the SB relations found for giants and super- 
giants are similar to the ones determined in the present 
paper for dwarfs and subgiants, especially their visible- 
infrared versions. This means qualitatively that any two 
stars of class I-V with similar magnitudes in two bands 
will present approximately the same angular diameters. 

8. Main Sequence stars as calibrators for 
long-baseline interferometry 

8.1. The need for small and nearby calibrators 

Interferometric observations are generally based on inter- 
leaved observations of a scientific target and a calibra- 
tor. The angular size of the calibrator is supposed to be 
known a priori, and the observed fringe contrast is used 
to estimate the instrumental transfer function (also called 
system visibility) . The catalogue of calibrators assembled 
by Cohen et al. (1999), and customized to interferometry 
by Borde et al. (|2002|l . consists mainly of K giants with 
angular diameters of about 2 mas. While this size is well 
adapted to short baseline observations (up to a few tens of 
meters in the infrared) , these stars are too large angularly 
to serve as calibrators for the hectometric baselines of the 
VLTI, the CHARA array (McAlister et al. l2~UUU|l or the 
NPOI (Armstrong et al.[1998). In addition, it is foreseen 
that shorter wavelengths will be implemented on the VLTI 
than the K band currently accessible with VINCI. For in- 
stance, the AMBER instrument (Petrov et al. [2000 ) will 
allow observations in the J band. The two-fold increase 
in angular resolution will naturally require significantly 



smaller calibrators than those in the Cohen et al. {1999) 
catalogue. 

A fundamental problem with distant stars is that the 
reddening corrections are uncertain. This means that it is 
highly desirable to use nearby stars as calibrators, located 
within a few tens of parsecs. In this respect, giant stars are 
not well suited due to their large linear dimensions, but 
dwarfs and subgiants are ideally suited to provide small 
and well-defined calibrators. 

Another advantage of Main Sequence stars is that 
their strong surface gravity results in a compact atmo- 
sphere and a well-defined photosphere. Their disk ap- 
pears sharper than that of the giants, for which the pre- 
cise definition of the limb darkened disk angular diameter 
at a level of less than 1 % can be difficult, in particular 
for the later spectral types. As an example, a discussion 
of the M4III giant -0Phe can be found in Wittkowski, 
Aufdenberg & Kervella l|2U04|l . 

8.2. Calibration precision vs. brightness 

It is possible to estimate the maximum angular size of 
calibrator stars in order to obtain a given relative pre- 
cision in the calibration of the interferometric efficiency. 
Fig. shows the achievable precision in the interferomet- 
ric efficiency using the (B, B — L) relation determined in 
Sect. 14.21 (rr < 1.0%), as a function of the angular diame- 
ter of the calibrator star, for four different baselines (100, 
200, 400 and 800m), in the H band. These baselines are 
representative of the existing or foreseen interferometers 
(Keck, PTI, VLTI, CHARA, NPOI and OHANA, sorted 
by increasing maximum baseline). The horizontal scale of 
Fig.Elcan be adapted for other wavelengths or other base- 
lines by scaling it linearly while maintaining constant the 
B/X ratio. 

If we now set a limit of 0.5% on the acceptable system- 
atic uncertainty in the interferometric efficiency, we can 
compute the apparent magnitude of the Main Sequence 
calibrators that should be used as a function of their color. 
The result is presented in Fig. das a function of the B — H 
color, for different baseline lengths and interferometric ob- 
servations in the H band. From this figure, it can be con- 
cluded that suitable calibrators for extremely long base- 
line observations will have to be faint. Let us consider 
the example of the OHANA interferometer (original idea 
proposed by Mariotti et al. ll99fij> . whose longest foreseen 
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0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 
Angular diameter (mas) 

Fig. 6. Precision achievable in the measurement of the 
interferometric efficiency as a function of the angular di- 
ameter of the calibrator, predicted using the (B,B — L) 
SB relation (< 1.0% dispersion). From left to right, the 
curves refer to baselines of 800, 400, 200 and 100 m, in the 
H band. 



baseline is 800 m. The H band magnitude of the calibra- 
tors necessary to obtain a relative systematic visibility er- 
ror of 0.5% will be between mn = 6 and 8, depending 
on the spectral type. This is rather faint, even for large 
aperture light collectors, but it is feasible with OHANA. 

As an alternative, it is possible to build (through time 
consuming observations) a secondary network of brighter 
and larger calibrators based on the small angular diameter 
ones, but there will always be a limitation attached to 
the fact that calibrators have to be observed in the first 
lobe of their visibility function. For OHANA, this sets a 
hard limit of ~0.5mas on the calibrator angular size, and 
even < 0.4 mas to obtain a visibility of at least 0.3. This 
corresponds to apparent magnitudes of ma = 5 to 7 in 
the H band, one magnitude brighter than the primary 
network. 

For the longest baseline of the VLTI (200 m), calibra- 
tor magnitudes between ran — 3 and 5 will be sufficient, 
clearly in the accessible domain of the AMBER beam 
combiner (Petrov et al. 12000(1 with the 1.8 m Auxiliary 
Telescopes (Koehler et al. 20<2J. The creation of a sec- 
ondary network of calibrators should therefore not be nec- 
essary. 

It should be stressed that the present conclusions re- 
garding magnitudes are not limited to dwarf stars, as gi- 
ants and supergiants follow comparable surface brightness 
relations. This means that the magnitude ranges defined 
above will be almost the same for other luminosity classes. 
A decisive advantage of dwarfs is that for the same appar- 
ent magnitude, they will be much closer than the more 
luminous classes, and therefore significantly less affected 
by interstellar extinction. 



9.0 



8.0 




1 3.0 



2.0 

l.o -I , , , , 1 

0.0 1.0 2.0 3.0 4.0 5.0 

B-H 

Fig. 7. Apparent magnitude in the H band of the calibra- 
tors suitable for obtaining a relative precision of 0.5% in 
the calibration of the interferometric efficiency, as a func- 
tion of the B — H color. From top to bottom, the curves 
refer to baselines of 800, 400, 200 and 100 m, in the H 
band. 



Table 8. Photometry (upper part) and predicted limb 
darkened angular diameters #ld (lower part) of the planet- 
hosting stars 51 Peg A and HD 209458 A. 





51 Peg A 


HD 209458 A 


mv 


5.50 ±0.01 


7.65 ±0.01 


Vtl K 


3.97 ±0.01 


6.31 ±0.03 


6 LD (K, V-K) 


0.689 ±0.011 mas 


0.228 ±0.004 mas 



8.3. Example of diameter prediction 

As a practical application, we have chosen the two stars 
51 Peg A (HD 217014) and HD 209458 A. The former hosts 
the first planet discovered around a solar type star (Mayor 
& Queloz 1995), and the latter presents planetary transits 
(Charbonneau et al. l2000fl . We selected these two stars be- 
cause they have been observed extensively using different 
techniques and did not show any large amplitude photo- 
metric variability. They therefore represent good examples 
of stable, well known stars, and are ideal candidates for 
the prediction of their angular size using the SB relations 
determined in the present paper. 

Tablc|Hlpresents the predicted angular diameters of the 
two stars for the (V, V—K) version of the SB relation. The 
my magnitudes are from Hipparcos (Perryman et al. H997[) 
for both stars, with an arbitrary error bar of ±0.01, while 
the K band infrared magnitudes were taken from Ducati 
et al. PKI2) for 51 Peg A, and from the 2MASS catalogue 
(Cutri et al. 12003) for HD 209458 A. 

For 51 Peg A, we obtain a predicted angular diame- 
ter of #ld = 0.689 ± 0.011 mas. The corresponding value 
for HD 209458 A is 6» LD = 0.228 ± 0.004 mas. These an- 
gular sizes can be translated into linear radii using the 
Hipparcos parallaxes (Perryman et al. I1997|l , 7751 p og a = 
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65.10±0.76 mas and 7thd209458A = 21.24 ±1.00 mas. We 
obtain i?5ip eg A = 1.138 ± 0.023 Rq and i?HD209458A = 
1.154 ± 0.059 R Q . HD 209458 A is a particularly interest- 
ing object, as Brown et al. I|2UU1|I have been able to esti- 
mate directly its linear radius through the deconvolution 
of the light curve of the transit. They obtain a value of 
^?hd 209458 a = 1.146 ± 0.050 Rq, in remarkable agreement 
with our (K, V — K) prediction. The bulk of the ±5% un- 
certainty comes from the error in the Hipparcos parallax, 
the relative error in the angular size being only ±2%. 

The direct measurement of the angular diameter of 
51 Peg A is within the capabilities of the existing very long 
baseline interferometers (several hundred meters), but this 
is not true for HD 209458 A. Its 0.228 mas size would re- 
quire baselines of more than 800 m to be resolved at visi- 
ble wavelengths (several kilometers in the infrared). Such 
baselines are not presently available or scheduled. And 
even so, the calibration of these observations would be ex- 
tremely difficult, as the calibrator would have to be very 
faint. More generally, carefully calibrated surface bright- 
ness relations are currently the only method to estimate 
precisely (±1%) the angular size of solar type stars fainter 
than my = 7. 

9. Conclusion 

The laws that we established between the angular size and 
broadband colors (or effective temperature) are strictly 
empirical. Our best relations present a very small intrin- 
sic dispersion, down to less than 1 %. They can be used 
to predict the angular sizes of A0-M2 dwarfs and A0-K0 
subgiants from simple, readily available broadband pho- 
tometry. On the one hand, Gray et al. H2003|l have re- 
cently published an extensive survey of the spectral prop- 
erties of nearby stars within 40 pc, including estimates of 
their effective temperatures. On the other hand, several 
large catalogues (2MASS, DENIS,...) provide high preci- 
sion magnitudes of these stars in the infrared. From the 
cross-comparison of these sources, the SB relations deter- 
mined in the present paper make it possible to assemble 
a catalogue of calibrators for interferometry that will be 
practically unaffected by interstellar extinction, multiplic- 
ity or circumstellar material biases. These resulting angu- 
lar diameter predictions will provide a reliable basis for the 
calibration of long-baseline interfcrometric observations. 
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